Cooperative Sparsity Pattern Recovery in Distributed 
Networks Via Distributed-OMP 



Thakshila Wimalajeewa and Pramod K. Varshney 
Department of Electrical Engineering and Computer Science, 
Syracuse University, Syracuse, NY 13244 
Email: twwewelw@syr.edu, varshney@syr.edu 



(N 

o 

(N 

> 
O 

00 
(N 



CZ2 



>: 

OV 



X 



Abstract — In this paper, we consider the problem of collaboratively 
estimating the sparsity pattern of a sparse signal with multiple measure- 
ment data in distributed networks. We assume that each node makes 
Compressive Sensing (CS) based measurements via random projections 
regarding the same sparse signal. We propose a distributed greedy 
algorithm based on Orthogonal Matching Pursuit (OMP), in which the 
sparse support is estimated iteratively while fusing indices estimated at 
distributed nodes. In the proposed distributed framework, each node has 
to perform less number of iterations of OMP compared to the sparsity 
index of the sparse signal. Thus, with each node having a very small 
number of compressive measurements, a significant performance gain 
in support recovery is achieved via the proposed collaborative scheme 
compared to the case where each node estimates the sparsity pattern 
independently and then fusion is performed to get a global estimate. 
We further extend the algorithm to estimate the sparsity pattern in a 
binary hypothesis testing framework, where the algorithm first detects 
the presence of a sparse signal collaborating among nodes with a fewer 
number of iterations of OMP and then increases the number of iterations 
to estimate the sparsity pattern only if the signal is detected. 

Keywords: Compressive sensing, Sparsity pattern recovery, 
multiple measurement vectors, distributed networks 

I. Introduction 

In the Compressive Sensing (CS) framework, a small collection 
of linear random projections of a sparse signal contains sufficient 
information for complete signal recovery (TJ — (3j . There is a con- 
siderable amount of work on the development of computationally 
efficient and tractable algorithms to recover sparse signals from CS 
based measurements obtained via random projections, for example 
in 1 4 1 — 1 12). However, there are several signal processing applications 
where complete signal recovery is not necessary. For example, in 
applications such as source localization in sensor networks [13 |, |14|, 
estimation of frequency band locations in cognitive radio networks 
1151 , localization of neural electrical activities from a huge number 
of potential locations in magnetoencephalography (MEG) and elec- 
troencephalography (EEC) for medical imaging applications |16|- 
1181 . sparse approximation 1 19 1, subset selection in linear regression 
|20], |21 1, and signal denoising [22], it is only necessary to estimate 
the locations of the significant coefficients of a sparse signal. 

A. Related work and our contribution 

The problem of sparsity pattern estimation of sparse signals has 
been addressed by several authors in the literature in the context of CS 
1 23 1-| 31 1. These studies focus on CS based sparsity pattern recovery 
with a single measurement vector. The achievable performance can 
be improved by having multiple observation vectors. Further, in 
distributed networks including sensor networks and cooperative cog- 
nitive radio networks, multiple measurements appear quite naturally 
since multiple nodes make observations regarding the same phe- 
nomenon of interest (Pol). Extensions of sparse signal recovery with 
multiple measurement vectors are considered in |32|, |33l assuming 
that the multiple measurements are available at a single location to 
perform the desired task. In such centralized settings, each node in 
the network has to transmit its observations along with the elements 
of the measurement matrix to solve the sparse signal recovery 
problem at the central fusion center. However, when performing 



inference tasks based on the observations collected at different nodes 
in distributed networks, the trade-off between the resource constraints 
(e.g. node power and communication bandwidth) and the achievable 
performance gain is a core issue to be addressed in many practical 
networks. Thus, distributed approaches for sparsity pattern recovery 
are desirable in many practical distributed networks where each node 
has limited computational and communication power/bandwidth. The 
problem of distributed sparsity pattern recovery is considered recently 
in the context of cognitive radio networks in 1341 , 1351 . There, 
decentralized consensus based algorithms for support recovery of 
sparse signals were proposed when each cognitive radio makes CS 
based measurements in cooperative cognitive radio networks. 

The work in this paper focuses on further reducing the computa- 
tional and communication burden at individual nodes in a distributed 
network while performing sparsity pattern recovery when each node 
in the network observes a noisy version of the same sparse signal. 
In contrast to transmitting raw observations along with measurement 
matrices to a central fusion center to perform centralized sparsity 
pattern recovery, we consider the case where each node tries to find 
an estimate of the sparsity pattern by collaboration and fusion. More 
specifically, we develop a greedy algorithm based on Orthogonal 
Matching Pursuit (OMP) where the indices of the sparse support 
are iteratively identified while fusing the estimated indices at each 
iteration. We show that, in the proposed distributed algorithm, each 
node has to perform less number of iterations compared to the sparsity 
index of the sparse signal to reliably estimate the sparsity pattern (in 
the centralized OMP algorithm, at least K iterations are required for 
sparsity pattern recovery where K is the sparsity index of the sparse 
signal). Moreover, in the proposed algorithm, each node transmits 
only one index at each iteration. Further, we extend the results and 
develop a joint algorithm to both detect the sparse signal and to 
perform the sparsity pattern recovery only if the sparse signal is 
detected. 

II. Problem Formulation and Background 

Consider a set of L distributed nodes making noisy measurements 
on a sparse signal. We assume that each sensor node obtains a 
M(< iV)-length measurement vector y; via CS based linear random 
projections. The measurement vector at i-th sensor node is given by, 



= Ais - 



(1) 



for / = 0, 1, • • • , L — 1. s is the N x 1 sparse signal and A; is 
the M x N random projection matrix at the i-th node. The noise 
vector V; at the l-th sensor node is assumed to be iid Gaussian with 
zero mean vector and the covariance matrix ct^Im where I m is the 
M x M identity matrix. 

When the signal s is sparse in the basis ^ such that s = where 
(3 contains only K « N number of significant coefficients, it has 
been shown that the randomized lower dimensional projections of 
the form l[T) can capture significant information of the sparse signal 
s. Assume that one has to estimate the sparsity pattern of a sparse 
signal which will give important information in many applications. 
For example, if the signal is sparse in Fourier domain, the locations 
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of non zero coefficients give the locations of significant harmonics 
which is important in spectrum sensing in the wideband regime. 

Let U be the support of the sparse coefficient vector f3 which is 
defined as, U := {i € {1, 2, • • • , N} \ (3{i) / 0} where f3(i) is the 
i-th element of f3 for i = 0, 1, • • • ,N—1. Then we have K = n(U) 
where n(S) denotes the cardinality of the set S. Further, let b be a 
iV-length vector which contains binary elements: i.e. 

= r i if (j® + o 

1 otherwise 

for i = 0, 1, • • • ,N — 1 and b be the estimated binary support vector. 

Sparsity pattern recovery with a single measurement vector can 
be performed via a CS reconstruction algorithm such as regularized 
least square algorithm (Lasso) |T), 1361 or OMP |8|. When multiple 
measurement vectors as in Q} are collected at a centralized location, 
the support of the sparse signal can be estimated, for example, 
by solving the simultaneous OMP algorithm as given in |32| 1 33 1 
in which the support can be directly estimated iteratively without 
reconstructing the complete sparse signal f3. 

III. Sparsity pattern recovery with multiple 

MEASUREMENT VECTORS: DISTRIBUTED ALGORITHM 
To implement a centralized sparsity pattern recovery algorithm 
based on the measurements collected at distributed nodes in a 
distributed network, it is required that each node transmits its M- 
length observation vector along with the elements of the measurement 
matrix to a central fusion center. Since transmitting all the information 
to a fusion center is not desirable in power and bandwidth constrained 
communication networks, we consider a distributed scheme in which, 
sparsity pattern estimation is performed via collaboration among 
nodes with less communication and with each node estimating the 
sparsity pattern. In 1341 , 1351 , several consensus based distributed 
schemes are proposed to estimate the support of sparse signals 
based on Lasso. These schemes estimate the full support set at each 
node and then perform fusion via several collaboration schemes. 
However, due to computational complexity, performing Lasso at 
power constrained sensor nodes may not be desirable. 

OMP is a greedy approach in which at each iteration, the location 
of one non zero coefficient or a column of G = A\& that participates 
in the measurement vector y is identified. More specifically, at each 
iteration, it picks the column of S which is most correlated with 
the remaining part of y. Then the selected column's contribution is 
subtracted from y and iterations on the residual are carried out. If we 
consider one sensor node, at each iteration, there are N — K possible 
incorrect indices that can be selected by the OMP algorithm. The 
probability of selecting an incorrect index at each iteration increases 
as the number of CS measurements at each node (A/) decreases. It 
has been shown that the OMP algorithm requires more measurements 
for signal reconstruction compared to optimization based (e.g. Basis 
Pursuit (BP)) algorithms. However, due to limitation of processing 
power at each node in many practical distributed networks, it is 
desired to have multiple nodes each processing a small number of 
measurements. Since all the nodes in the network observe the same 
sparse signal, it is highly likely that the estimate of the index at 
multiple nodes is the same at a given iteration of OMP. To reduce 
the error in incorrectly selecting an index at each iteration of OMP 
with less number of compressive measurements, we propose to fuse 
the indices estimated by each node during a given iteration by 
collaboration among the nodes in the network. It is worth mentioning 
that, in the proposed approach, a node in the network may find a set 
of indices (instead of one index as in the conventional OMP) that 
correspond to non zero coefficients via collaboration, especially when 
the number of distributed nodes is closer to or greater than the sparsity 
index. Thus the proposed algorithm can be terminated with less 
number of iterations compared to K at each node. It should be further 
noted that, in the proposed OMP based algorithm, at each iteration, 



only one index needs to be transmitted by each node for fusion, thus 
reducing communication cost compared to that with distributed Lasso 
versions proposed in [34], [35|. We dub the proposed algorithm as 
'Distributed and collaborative OMP (DC-OMP)'. 

Distributed and collaborative OMP for sparse support estimation: 
Define Mi to be the set containing the neighboring nodes of the 
Z-th node (including itself). As defined in Subsection [TTJ let U be the 
support set which contains the indices of non zero coefficients of the 
sparse signal and lAi be the estimated support set at the Z-th node. 
Further, let Qi — Ai^ and 8i t i be the i-th column of the matrix 
Oi. Oi(A) denotes the submatrix of O; which has columns of 0; 
corresponding to the elements in the set A for A C {0, 1, • • • , N — 
1}. |.| denotes the absolute value while denotes the l p norm. 
Further, n(S) denotes the cardinality of the set S as defined before. 

In the proposed DC-OMP algorithm which is stated in Algorithm 
[T] once the Z-th node finds an index A;(t) (which corresponds to 
the column that is most correlated with the remaining part of yi) 
at iteration t by performing step 2 in Algorithm Q] it is exchanged 
among the neighborhood Mi. Subsequently, each node will have the 
index that the nodes in its neighborhood obtained from step 2. By 
fusion, each node selects a set of indices (from n(Mi) number of 
indices) such that most of the nodes agree upon (more details of this 
fusion are provided in Subsection llll-lt . Note that, in this step several 
indices can be selected and thus, the number of iterations required to 
estimate the support that each node has to perform can be less than 
the sparsity index K. 



Algorithm 1 Distributed OMP for sparsity pattern estimation 
At Z-th node: 

1) Initialize t = 1, W;(0) = 0, residual vector r;,o = y; 

2) Find the index A; (t) such that 

\ M I < *l,t-i,0l,u > I 
Xi(t) = arg max -f -y 

u=i,— ,n ||ri,t-i||a 

3) Update the index set A ; *(t) via local communication: X*(t) = 
fi(Xi(t),{Xi(t)},i £ Mi), as discussed in subsection IIII- II 

4) Set Ui(t) = Ui(t - 1) U {\l{t)} and l t = n(U,(t)) 

5) Compute the projection operator Pi(i) = 

@i(Ui(t)) (ei(Wj(t)) T ei(Wi(*))) _1 e,(w,(t)) T . update 

the residual vector: ri.t = (I — P;(t))y; 

6) Increment t = t + 1 and go to step 2 if It < K, otherwise, set 



1) Performing step 3 in Algorithm [7} To perform step 3 in 
Algorithm Q] we propose the following procedure. 
Case P. Consider the case where the Z-th node can broadcast its 
estimated index at each iteration to the rest of the nodes in the 
network, i.e. Mi — M where M contains the indices of all 
the nodes in the network. This is a reasonable assumption when 
there are only a small number of nodes in the network (e.g. 
cognitive radio networks with a 5 — 10 cognitive radios). During 
each iteration t of the distributed OMP algorithm, the Z-th node 
broadcasts A; (£). Consequently, the Z-th node receives the estimates 
\i(t)'s for i G M from the whole network. Further, let c{t) 
be a L-lengfh vector that contains all the indices estimated at L 
nodes from step 2 during the t-th iteration (i.e. values of Xi(t) for 
i = 0, 1, • • • , L — 1). At t-th iteration, X*(t) is updated as follows: 
Check whether there are indices with more than one occurrences 
(i.e. whether there is any index in the vector c(t) that occurs more 
than once). If yes, such indices are put in the set \1(t) (such that 
A*(t) = {set of indices which occur more than once}. If no, 
(i.e. there is no index obtained from step 2 such that any two or more 
nodes agree with, so that all L indices in c(t) are distinct), then select 
one index from c(t) randomly and put in X*(t). In this case, to avoid 
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the same index being selected at subsequent iterations, we force all 
nodes to use the same index. 

It is noted that when there are approximately equal or more 
distributed nodes compared to sparsity index K, the vector c(t) has 
at least one set of two indices with the same value, thus, X*(t) is 
updated appropriately most of the time. Further, in this case, since 
each node has the indices received from all the other nodes in the 
network, every node has the same estimate for U after algorithm 
terminates. 

Case II: Next, we consider the case where Mi C M; i.e. each 
node communicates its estimated index in its neighborhood which has 
fewer number of nodes compared to all the nodes in the network. 
Then, similar to the above case A*(t) is found based on c;(t) as 
the indices which have the most occurrences from ci (t) where a (t) 
contains the indices received by the l-lh node from its neighborhood 
at t-lh iteration. However, in this case, since Z-th node does not receive 
the estimated indices from the whole network, at each iteration, 
different nodes may agree upon different indices (however, we see 
experimentally that for a relatively large size of Mi, all the nodes 
have the same estimated index set at the end). When two neighboring 
nodes agree upon two different indices at t-lh iterations, there is a 
possibility that one node selects the same index at a later iteration 
than t. To avoid the l-lh node selecting the same index twice, we 
perform an additional step; i.e., to check whether X*(t) determined 
as above is in Ui{t - 1). If X*{t) € Ui(t - 1), set A*(t) = Ai(t), 
otherwise update A* (t) similar to the procedure described in Case I. 

IV. Sparse signal detection and sparsity pattern 

ESTIMATION 

Next, we consider the case where one has to detect whether the 
sparse signal is present and estimate the sparsity pattern only if the 
signal is present. We consider the following binary hypothesis testing 
problem, 



Hi : yi = A ; s + V;; 
Ho ■ y; = vj. 



(2) 



for I = 0, 1, • • • , L — 1 and where under hypothesis Hi the sparse 
signal is present. In the following, we extend the collaborative 
algorithm discussed above to first detect the sparse signal and then to 
estimate the sparsity pattern without ever reconstructing the signal. 
Further, we assume that L > K. 

The idea is based on the properties of the OMP algorithm. If the 
signal is not present in the model l[2}, it is very unlikely that two nodes 
in the network will select the same index of the support set at any 
given iteration based on the step 2 in the OMP algorithm presented 
in Algorithm [T] However, when the signal is present (i.e. hypothesis 
Hi), the probability that two nodes select the same index at each 
iteration is higher as the number of nodes is close to or greater than 
the sparsity index K. That is because, when the signal is present and 
in relatively high signal-to-noise ratio region, the column index of the 
projection matrix Qi which is most correlated with the remaining part 
of the observation can be estimated at the Z-th node as one of the 
index of the true sparse support, and multiple nodes can get the same 
index during a given iteration. Taking this information into account, 
we extend the algorithm such that it first detects the sparse signal with 
fewer number of iterations and increase the number of iterations to 
find the sparsity pattern only if the signal is detected as described in 
Algorithm |2] 

2) Updating ii n dex in step 4: Step 4 in Algorithm [2] is performed 
as given below. At t-th iteration, q (t) contains all the indices received 
by the Z-th node from its neighborhood. The function unique(ci(t)) 
gives the number of distinct indices of the support set in c(t). If 
all the indices in c;(t) are different from each other, unique(ci(t)) 
equals to the number of nodes in the neighborhood of the l-lh node 
including itself. If there are any two indices in ci(t) with the same 



Algorithm 2 Distributed OMP for sparse signal detection and sparsity 
pattern estimation 
At l-th node: 

1) Initialize t = 1, W;(0) = 0, residual vector r;,o = yi, iindex = 


2) Find the index A; (t) such that 

' < rj (t _i,0i lU > | 



arg max 

(*>=!,••• ,N 



\ri,t-i 2 



3) Update the estimated index set A; (t) via local communication: 
A*(t) = fi(Xi(t), {Xi(t)},i £ Mi), as discussed in subsection 

4) Update i index : 

• if unique(ci(t)) = n(Mi), iindex = iindex + 

• else ( unique(ci(t)) < n(Mi)), i in dex = iindex + 
P(ci(t)) 

as discussed in Subsection IIV-21 

5) Perform signal detection decision when t — Ko 

• If t = A'o and ii n dex(t) > lo, decide Hi and go to step 
6. Avoid steps 4 and 5 in subsequent iterations 

. If t = K and i in dex(t) < I decide H , set Ui(t) = 
and go to step 9 

6) Set Ui(t) = Ui(t - 1) U {Xi{t)}, and l t = n{Ui{t)) 

7) Compute the projection operator Pi(t) = 
&i(Ui{t)) (Oi(tii(t)) T Gi{tii{t))y @i(Ui(t)f. 
the residual vector: r;. t = (I — P;(t))yj 

8) Increment t = t + 1 and go to step 2 if l t < K, 

9) set Ui =Ui{t) 



Update 
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Fig. 1. Performance of the sparsity pattern recovery with distributed OMP 

algorithm: (left) the probability of correctly recovering the sparse support 

(Po = Pr(h = b)) vs M/N (right) the percentage of the support correctly 
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value, we set the value of p(ci(t)) as the number of such indices. 
After performing Ko (which is less than K) number of iterations, 
if iindex hi step 4 in Algorithm [2] is very small (less than Io where 
lo *C K), the algorithm decides that no sparse signal is present and 
terminate the process resulting in the null set as the estimated support 
set. If ii nd ex > lo, it decides that the sparse signal is present and 
continues estimating the support set similar to Algorithm Q] 

V. Simulation Results 

In this section, we present some simulation results to demonstrate 
the performance of proposed algorithms for sparsity pattern recovery 
based on distributed CS based measurements. In the following, 
we assume that the entries of each projection matrix A; for I — 
0, 1, • • • , L — 1 are drawn from a Gaussian ensemble with mean zero 
and variance i, 

To compare the performance of the proposed Algorithm Q] with 
other approaches, we consider two benchmark cases, (i). Distributed 
OMP with no collaboration: in this case, each node performs OMP 
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K=10, L=10, N=256, SNR=1 7.3227 dB 
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Fig. 2. Number of iterations of OMP at each node with the proposed Algo- 
rithmQJvs M/N; N = 256, K = 10, L = 10, 7 = = 17.3227(iB, 



independently to obtain the support set estimate Ui(t), i.e., the step 
3 is eliminated in Algorithm QJ and set \*{t) = \i(t). To fuse 
the support sets estimated at individual nodes, each node transmits 
its estimated support to a fusion center and performs a majority 
rule based fusion scheme to obtain a global estimate of U. (ii). 
Simultaneous OMP (S-OMP) |32|: S-OMP algorithm is carried out 
using all the raw observations at the fusion center. Thus, this scheme 
requires each node to transmit its M -length observation vector as 
well as the projection matrix A; to the fusion center. 

In Figures QJ and [2] we illustrate the performance of the sparsity 
pattern estimation based on proposed DC-OMP as in Algorithm [Tj 
in terms of different performance metrics. In both figures, we let 

K = 10, N = 256, L = 10, and SNR at each node, defined as 

11 11 2 

7 = ILlij = 17.3227dP. Also, in Figures QJ and [2] we assume 
that Mi = M as considered in case I in Subsection IIII-ll Then 
the estimated support set at each node based on Algorithm [Tj is 
the same. In Fig. [TJ by performing 10 4 runs and averaging over 
10 trials, we plot the probability of correctly recovering the full 
support set, Pd = Pr(b = b) (left) and the percentage of 
the support set that is estimated correctly (right) vs M /N where 
M is the number of compressive measurements at each node. It 
can be seen from Fig. QJ that, at relatively small values of M/N, 
the proposed algorithm outperforms D-OMP with no collaboration. 
In resource constrained distributed networks, especially in sensor 
networks, it is desirable to perform the desired task by employing 
less measurement data (i.e. with small M) at each node distributing 
the computational complexity among nodes to save the overall 
node power. Thus, fusing the estimated indices of the non-zero 
coefficients at each iteration of the OMP algorithm ensures a higher 
performance in exact sparsity pattern recovery compared to that when 
OMP is performed at each node independently. However, as M/N 
increases, Pr(b = b) of both algorithms converge to 1, since when 
the number of compressive measurement at each node increases, 
OMP (with or without collaboration) works better and recovers the 
sparsity pattern almost exactly with a single measurement vector. 
It has been shown in (8j, (37 1 that OMP requires approximately 
M > (1 + e^.K' log TV, e > measurements for reliable sparsity 
pattern recovery in the noise free case. Thus, even with a very small 
M at each node, having L number of nodes, S-OMP achieves this 
limit and provides a significant performance gain compared to the 
proposed algorithm at very small M/N values. However, S-OMP 
requires a considerable communication overhead compared to the 
proposed algorithm. Further, in the proposed algorithm, each node 
has the same estimator at the end in contrast to the centralized S- 
OMP. In Fig. [TJ we further plot the percentage of support that is 
correctly recovered, for example, at ^ w 0.1, the proposed algorithm 
correctly recovers approximately 75% of the support while D-OMP 
with no collaboration recovers only about 30% of the support. Thus, 
from both sub figures of Fig. [TJ significant performance gain is 







✓* 












/ 








































' -"Sparse signal detection (algo 2) 
— Soars ty patten dr:oct on (a go 2) 
* Soars iy patrcyn dr:oct on (a go 1) 





0.06 
EO-075 




















rt 0.07 
03 

W0.065 










ni 

III 0.06 










o 

>G.055 










•••sparsity pattern 


-QO.05 






• sparse signal 


_Q 










00.045 

Q_ 

0.04 






•..•••»-... 









1 0. 


5 


2 0. 



M/N 



M/N 



Fig. 3. Performance of the sparse signal detection and sparsity pattern 

- = Ml = 

17.3227c(B; (left) Probability of detection, (right) probability of false alarms 



recovery with Algorithm [2] N = 256, K = 10, L 



10, 



observed via the proposed algorithm compared to D-OMP with no 
collaboration with small M, which is the more desirable scenario 
in resource constrained distributed networks. To further illustrate the 
efficiency of the proposed algorithm, in Fig. [2] we plot the number of 
iterations of the DC-OMP algorithm that each node has to perform 
in recovering the sparsity pattern. It is observed from Fig. [2] that as 
M /N increases, the proposed algorithm estimates the sparsity pattern 
reliably by executing only ~ K/2 number of iterations at each node. 
When M/N increases, as observed from Fig. [TJ the performance of 
both DC-OMP and D-OMP with no collaboration converges to the 
same level but DC-OMP requires very small number of iterations at 
each node to achieve that performance compared to that with D-OMP 
with no collaboration which requires K number of iterations at each 
node irrespective of the value of M/N . 

In Fig. [3] we illustrate the performance of Algorithm[2]for detecting 
the sparse signal before estimating the sparsity pattern. We plot 
the performance of sparse signal detection as well as the sparsity 
pattern estimation in Fig. [3] For sparse signal detection, probability 
of detection and the false alarm are given by Pf> = Prob(S = 
and Pp — Prob(S = l\Ho), respectively where 8 is the detection 
decision. For sparsity pattern detection, the probability of detection 
is given by P% = Pro6(b = b 1 \Hi)Prob(Hi) + Prob(b = 
h°\Ho)Prob{T-La), where we redefine the variables such that b 1 
is the binary support of the signal s (i.e. the support under Hi) 
while b° denotes the vector with all zeros (binary support under T-Lq). 
Similarly the probability of false alarm, is given by Pp = Prob(h = 
b 1 |Ho)Pro6(Hi) + Pro6(b = b°\Hi)Prob(Ho). Further, in Fig. [3] 
we use the same values for the parameters TV, K, L and 7 as used in 
Figures QJ and [2] and Ko = 3 and 7o = 2. From Fig. [3] it is seen that 
Algorithm [2] reliably detects the sparse signal even with a very small 
value of M/N, and the performance of the sparsity pattern recovery 
after detecting the signal has performance that is close to that when 
the sparsity pattern recovery is done as in Algorithm QJ (where it is 
known a priori that the signal is present). 

VI. CONCLUSION 

In this paper, we addressed the problem of recovering a com- 
mon sparsity pattern based on CS measurement vectors collected 
at distributed nodes in a distributed network. A distributed greedy 
algorithm based on OMP is proposed to estimate the sparsity pattern 
via collaboration in which each distributed node is required to 
perform less number of iterations of the greedy algorithm compared 
to the sparsity index. When it is not known a priori that a sparse 
signal is present or not, the algorithm was extended to perform 
detection of the sparse signal with a fewer number of iterations before 
completely recovering the sparsity pattern. The proposed algorithm 
is shown to have significant performance gains compared to that 
with each node performing OMP independently and then fusing the 
estimated supports to achieve a global estimate. Complete theoretical 
analysis of the algorithm will be considered in a future work. 
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